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Abstract. We study the relativistically expanding 
electron-positron pair plasma formed by the process of 
vacuum polarization around an electromagnetic black hole 
(EMBH). Such processes can occur for EMBH's with mass 
all the way up to 6 • lO^Af0. Beginning with a ideal- 
ized model of a Reissner-Nordstrom EMBH with charge 
to mass ratio ^ = 0.1, numerical hydrodynamic calcu- 
lations are made to model the expansion of the pair- 
electromagnetic pulse (PEM pulse) to the point that the 
system is transparent to photons. Three idealized special 
relativistic models have been compared and contrasted 
with the results of the numerically integrated general rel- 
ativistic hydrodynamic equations. One of the three mod- 
els has been validated: a PEM pulse of constant thick- 
ness in the laboratory frame is shown to be in excellent 
agreement with results of the general relativistic hydro- 
dynamic code. It is remarkable that this precise model, 
starting from the fundamental parameters of the EMBH, 
leads uniquely to the explicit evaluation of the parameters 
of the PEM pulse, including the energy spectrum and the 
astrophysically unprecedented large Lorentz factors (up 
to 6 • 10^ for a lO^M© EMBH). The observed photon en- 
ergy at the peak of the photon spectrum at the moment 
of photon decoupling is shown to range from 0.1 MeV to 
4 MeV as a function of the EMBH mass. Correspondingly 
the total energy in photons is in the range of 10^^ to lO'^^ 
ergs, consistent with observed gamma-ray bursts. In these 
computations we neglect the presence of baryonic matter 
which will be the subject of forthcoming publications. 
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1. Introduction 

At present there is a need for viable source models capa- 



ble of producing the large energy 
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10^"* ergs) 



necessary to po wer distant gamma-ra y bursts such as 



the production of such large bursts based on results of cal- 
culations of the energy emission from vacuum polarization 
processes around an EMBH. 

Damour & Rufhni (1975) indicated a very precise 
mechanism to explain the origin of gamma ray bursts: the 
vacuum polarization of a Kerr-Newman EMBH via the 
Heisenberg-Euler-Schwinger process. They showed how 
these processes, occurring for EMBH's all the way up to 
6 • IO^Mq, can radiate a large fraction of the extractable 
energy of the EMBH, in t he sense of the Black Holes re - 
versible transformations (IChristodoulou fc Rufhni 19711) . 
Recently, after the discovery of the gamma ray afterglow 
by t he Beppo Sax satellite and their optical identification 
(e.g. iFrontera fc Piro 199S|) , this model has been reconsid- 



ered by Rufhni (1998) who introduced the new concept of 
the " dyadosphere" of an EMBH and by Preparata et al. 
(1998a, 1998b) who have given details of the spatial density 
and energy distribution in the electron-positron plasma 
created in the dyadosphere. The present paper examines 
the time evolution of such an electron-positron plasma by 
assuming an EMBH in vacuum and calculates the conse- 
quences. 

Past work on black hole formation by Wilson (1975, 
1977) has shown that highly charged EMBH's with charge 
to mass ratio ^ = Q/{mVG) « 0.1, surrounded by an 
oppositely charged magnetosphere, may form during the 
collapse of rotating magnetized stars. For a low mass star 
to collapse and form a highly charged black hole, sur- 
rounded by an oppositely charged magnetosphere, strong 
initial fields and rapid rotation are required. The theoreti- 
cal background for such an analysis has been presented in 
Rufhni & Wilson (1975), and Damour et al. (1978). The 
rece nt discovery of the soft g a mma-ray repeater SGRl 806- 
20 (iKouveliotou et al. 1998i [Baring fc Harding 199^) in- 



GRB971214 (e.g. [Frontera fc Piro 1999|) and GRB990123 



( Bloom et al. 19991 ). In this paper we present a model for 



dicates that very high magnetic fields may be present in 
precollapse stars of a few solar masses. 

The calculations by Wilson assumed space was 
filled with a perfect conductor so that the magneto- 
hydrodynamic model could be used. In the region above a 
newly formed EMBH the density is falling rapidly. Even- 



2 



Remo RufBni et al.: On the Pair Electromagnetic Pulse of a Black Hole with Electromagnetic Structure 



tually, due to the paucity of charge carriers, the assump- 
tion of the plasma as a perfect conductor must fail, pos- 
sibly leading to regions where the electric fields will be 
sufficiently high so that pair creation will ensue. Similar 
scenarios can be envisaged for gravitational collapse of 
much larger masses M ~ 10* — IQ^Mq occurring in galac- 
tic nuclei. In the present paper we begin with the ideal- 
ized situation of an isolated, Reissner-Nordstrom EMBH 
with charge to mass ratio ^ = 0.1. The case of a Ker 



Newmann EMBH is presented in pamour fc Ruffini 1999 . 
The vacuum polarization processes are computed, follow- 
ing the approach outlined in Rufhni (1998), Preparata et 
al. (1998), Ruffini (1999a), Ruffini et al. (1999), Ruffini 
(1999b), and they give the specific initial data of the dya- 
dosphere. The evolution of the PEM pulse formed by the 
vacuum polarization processes is followed by a numerical 
calculation to the point where the system is transparent to 
photons. We then analyze the final gamma-ray signal that 
would be observed from this idealized model of an EMBH 
in vacuum. Extensions of this model such as the expansion 
of the PEM pulse in the presence of the remnant baryonic 
matter left over in the collapse process will be the subject 
of a forthcoming publication Ruffini et al. 1999b. Some 



consequences for the interpretation of the observed burst 



and afterglow are presented in Preparata & Ruffini 1999 



and Preparata, Ruffini & Xue 1999 



2. Formation of e~^e pairs from an EMBH 

We recall the relevant formulas from Ruffini(1998) and 
Preparata et al. (1998a, 1998b): 

The Christodoulou-Ruffini mass formula 



for a Reissner-Nordstrom 
( IChristodoulou fc Ruffini 1971] ) 



E = Mc^ 



2r+-- 



fi2 

S = Anrl = IGtt ( 



with 



Vgq 

r+c^ 



< 1, 



black hole gives 

(1) 
(2) 

(3) 



where Mjr is the irreducible mass, r+ = 2Mir is the hori- 
zon radius, S is the horizon surface area, and maximally 
charged (Qmax — 2^/GMir) black holes satisfy the equal- 
ity in Eq.(||). 

The dyadosphere is defined as the region outside the 
horizon of a EMBH where the electric field exceeds the 
critical value for spontaneous e+e^ pair creation, £c = 



( Peisenberg fc Euler 193l| , pchwinger 1951D , where 
m and e are the mass and charge of the electron. By in- 
troducing the dimensionless mass and charge parameters 
fi = ^ — < 1 for Reissner-Nordstrom black 



holes, the horizon radius may be expressed as 



GM 



r+ = 



GAP 



= 1.47 X 10V(1 + Vl -^^) cm. 



(4) 



The outer limit of the dyadosphere is defined as the radius 
Tds at which the electric field produced by the EMBH 
equals the critical field. The radius of the dyadosphere can 
be expressed in terms of the Planck charge qp = {fic)^ and 
the Planck mass m-p — (^)^ in the form 



mcj \ P J \m) \qpj \y^M J 



1.12 



V c2 y 

IOVmC cm. 



(5) 



which well illustrates the hybrid gravitational and quan- 
tum nature of this quantity. The radial interval r+ < r < 
Tds describes the region in which e^e~ pairs are produced 



(see Fig.[l] of Preparata, Ruffini & Xue 1998a). 

In a very short time 0{:;^)~ a very large number 
of pairs is created in the dyadosphere. The local number 
density of pairs created is computed as a function of radius 



ne+e-{r) 



Q 



A'Kr'^ 



1 



and the energy density is given by 
.(r).^fl ' ^ 



87rr4 



(6) 



(7) 



(see Figs. [2] & [3] of [Preparata, Ruffini fc Xue 1998aD . 
The total energy of pairs converted from the static electric 
energy and deposited within the dyadosphere is then 



2 r_ 



-(1 



rds 



1+ 
rds 



(8) 



In the limit of — 



0, Eq.d) leads to Elf -> i^, 
which coincides with the energy extractable from EMBH's 
by reversible processes {Mir = const.), namely E — Mir = 
i2!(see Fig. [4] of [Preparata, Ruffini fc Xue 1998a| ). 

Due to the very large pair density given by Eq. (^ and 
to the sizes of the cross-sections for the process e^e^ <-> 
7 -|- 7, the system is expected to thermalize to a plasma 
configuration for which 



N, 



7 ~ ^'P 

and reach an average temperature 



TT'tot 

3Np ■ 2.7' 



(9) 



(10) 



where iV„ 



is the number of e+e -pairs created 
in the dyadosphere and k is Boltzmann's constant (see 
Fig. [4] of Preparata et al. 1998b). 
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1000 ^ 




100 1000 
Mass/Mg 

Fig. 1. The total energy in e+e^ pairs deposited within 
the dyadosphere of a black hole with mass M and charge 
Q = O.lQmax equal to 10% of maximal black hole charge. 




100 1000 
Mass/Mg 

Fig. 2. An average energy per pair versus mass is given 
for a black hole with mass M and charge equal to 10% of 
maximal black hole charge: ^ = 0.1. 



In conclusion, we stress that (i) the locality of e+e^ 
pair creation is an important feature of the dyadosphere, 
and (ii) e~^e~ pair creation occurs only in a finite range 
of the parameters /i and ^ and a unique spectrum of the 
created pairs can be calculated in terms of these two pa- 
rameters /u,^ (iii) Eqs.(|l]-||) show that the pair-creation 
process is capable of extracting large amounts of energy 
from an EMBH with an extremely high efhciency (close 
to 100%). 

In Wilson (1975, 1977) a black hole charge of the order 
Q = 0.1(5max(C — 0-1) was formed. Thus, we henceforth 
assume an EMBH with ^ = 0.1 for our detailed numerical 
calculations. In Fig.|^ a plot of the total energy Eq.(||) 
versus mass is given for the case ^ = 0.1 and in Fig. a 
plot of the average energy per pair (relating to Eq.(O)) 
versus mass is given also for ^ = 0.1. 



3. General relativistic hydrodynamic equations 

In order to model the radial evolution of the e+e^-pair 
and photon plasma fluid created in the dyadosphere of an 
EMBH, we need to discuss the relativistic hydrodynamic 
equations governing such an evolution and the associated 
rate equation for the pairs. 

The metric for a Reissner-Nordstrom black hole is 



2 • 2 

r sm 



2 _ 2GM 



(11) 
and 



ds^ = -gtt{r)df + grr{r)dr'^ + r'^dO'^ + 

where gtt{r) = 

grr{r) = a{ry^. 

The plasma fluid of e"''e~-pairs, photons and baryon 
matter in the curved spacetime given by Eq.(pT|) is de- 



scribed by the covariant stress-energy tensor given by 
( [Weinberg 197^ ) 

T'^" = pg^^^ + {p + p)U''U'' + AT^", (12) 

where p and p are respectively the total proper energy 
density and pressure in the comoving frame. C/^ is the 
four-velocity of the plasma fluid. AT'^^ accounts for dissi- 
pative effects due to heat conduction, shear and bulk vis- 
cosity. In this paper, as a flrst approximation, we assume 
the plasma fluid of e'^e~ pairs, photons and baryon matter 
to be a simple perfect fluid in curved spacetime and ne- 
glect the dissipative terms AT'^'^ in Eq. (p^) . In spherical 
symmetry we have 



9tt{uy + grriuy 



1 



(13) 



where and [/* are the radial and temporal contravari- 
ant components of the 4- velocity. The conservation law for 
the baryon number can be expressed as a function of the 
proper baryon number density rt^ as follows 



1 



(14) 



The energy-momentum conservation law of the plasma 
fluid is 



1 dg^x 

2 dxt" 



9 



ip- 







(15) 



The radial component reduces to 



dt 



dr 



p)U''Ur 
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dgtt 
dr 



t\2 



dgr 



r\2 



dr 







(16) 



The component of the energy-momentum conservation 
equation (nsh along a flow line is 



1 , 1 



g-H9^pun,.-P9-H9~'un, 

1 



p 



= . 



(17) 



Defining the total proper internal energy density e and the 
baryon mass density pB in the comoving frame 



e^p- Pb, 



Pb = nBrriBC 



(18) 



and using the law of baryon- number conservation (p^, 
from Eq.(n^) we have 



{eUn;.+pK 







(19) 



Recalling that ^ = VU.^^, where V is the comoving vol- 
ume and T is the proper time for the fluid, we have along 
each flow line 



dr 



dV 



■P- 



dr dr 



dE dV 

P— = 



dr 



(20) 



where E = Ve is total proper internal energy of the plasma 
fluid. 

Finally, we represent the equation of state by a thermal 
index T{p, T) 



1 



P 
e 



(21) 



We now turn to the analysis of e+e~ pairs initially 
created in the dyadosphere. Let rig- and ng+ be the proper 
densities of electrons and positrons. We have 



(22) 



The rate equation for electrons and positrons gives 
= CT!7[ne-(T)ne+(r) - ne-ne+] , 



(23) 



where a is the mean pair annihilation-creation cross- 
section, V is the thermal velocity of e^, Wv is the average 
value of av and rigi (T) are the proper number-densities 
of electrons and positrons that thermalize with photons 
through the process e+e^ ^ 7 + 7, ^e±(r) = n-y{T) for 
T > iJieC^, where n^(T) is the number-density of pho- 
tons. The proper number-densities rigi (T) are given by 
appropriate Fermi integrals with zero chemical potential, 
at the equilibrium temperature T . The equilibrium tem- 
perature T is determined by the thermalization processes 



occurring in the PEM-pulse with a total proper energy- 
density /9, whose evolution is governed by the hydrody- 
namic Eqs.(|l|,|l6|,0). We have 



P = P-t + Pe+ + Pe- + 



P\- 



PB, 



(24) 



where p-y = aT^ is the photon energy-density and 
is the energy-density of the electrons stemming from the 
ionization of baryonic matter, which is given by a Fermi 
integral with an appropriate chemical potential p^ at the 
equilibrium temperature T, and pB — mBC^riB, where ub 
is the baryon number-density, since baryons and ions are 
expected to be non-relativistic in the range of temperature 
T under consideration. p^± is the proper energy-density of 
electrons and positrons of pairs deflned by 



Pe± 



Pe± {T), 



(25) 



where n^± is obtained by integration of Eq.(^3|) 
and pg± (T) is the proper energy-density of electrons 
(positrons) obtained from zero chemical potential Fermi 
integrals at the equilibrium temperature T. Having im- 
plicitly deflned the equilibrium temperature T by Eq. ( p^ , 
we can also analogously evaluate the total pressure 



P = Pj+Pe+ +Pe- 



pI 



PB, 



where p-y is the photon pressure, Pe± defined by 

n„± 



Pe± 



n^± (T) 



Pe±{T), 



(26) 



(27) 



where the pressures p^± (T) are determined by zero chem- 
ical potential Fermi-integrals at the equilibrium temper- 
ature T, and is the pressure of the ionized electrons, 
evaluated by an appropriate Fermi integral of non-zero 
chemical potential p^. In Eq.(|2^), the ion pressure is 
negligible by comparison with the pressures p^ ^± (T) . 
Finally, using Eqs.( |2^ , |2^ ), we compute the thermal index 
F of the equation of state (|2l|). 

It is clear that in order to be correct, the entire set of 
equations considered above, namely Eqs. (|l^ , [l^ , |l7|) with 
equation of state given by Eq.(^) and the rate equa- 
tion ( |2^ ) have to be integrated fulfilling the total en- 
ergy conservation for the system. The boundary condi- 
tions adopted here are simply purely ingoing conditions 
at the horizon and outgoing at radial infinity. 

The calculation is initiated by depositing a proper en- 
ergy density (0) between the Reissner-Nordstrom horizon 
radius r_|_ and the dyadosphere radius r^s- The total en- 
ergy deposited is given by Eq.(^). The calculation is con- 
tinued as the plasma fluid expands, cools and the e+e^ 
pairs recombine, until it becomes optically thin: 



dr{n^± + ZnB)<JT — 0(1), 



(28) 



where ctt — 0.665 x 10 ^"^cm^ is the Thomson cross- 
section is , Z is the average atomic number of baryon 
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matter and the integration is over the radial interval of 
the PEM-pulse in the comoving frame. At this point the 
energy is virtually entirely in the form of free-streaming 
photons and the calculation is stopped. In the following 
computations we shall assume that baryonic matter is ab- 
sent. The case of presence of baryonic matter, left over in 
a remnant during the process of gravitational collapse, is 



discussed in Ruffini et al. 1999b 



4. Numerical integration of general relativistic 
hydrodynamic equations 

We use a computer code (Wilson et al. 1997, 1998) to 
evolve the spherically symmetric hydrodynamic equations 
for the baryons, e'*'e~-pairs and photons deposited in the 
dyadosphere. 

We define the generalized Lorentz factor 7 and the 
radial coordinate velocity 



7 = v/1 + U-Ur 



From Eqs.(|ll|, |l 
1 

9tt 



we then have 
1 



(30) 



Following Eq.(|l8|), we also define 

E = ej, D = pBj, and p = 

so that the conservation law of baryon number ( p^ can 
then be written as 

dD ct d r^ 
— = - — — {—DV). 
dt r^ dr a 

Consequently, Eq.(^^ acquires the form. 



(31) 



(32) 



dE 
'dt 



a d r^ 
^^C-EV^ 
r^ or a 



P 



'dt 



a d r^ 



dr 



a 



(33) 



Defining the radial coordinate momentum density 

Sr = a{p + p)U^Ur = iD + rE)Ur, (34) 

we can express the radial component of the energy- 
momentum conservation law given in Eq.(p^) by 

dSr 
dt 



r-^ or a or 



a 



ip + p) 



tt 



dg_ 

dr 



t\2 



dgr 



dr 



r^ or a or 



M 



D + TE 

7 



(35) 



In order to determine the number-density of e'^e~ 
pairs, we turn to Eq.(p3|). Defining the e"'"e~-pair coordi- 
nate density A'^g± = ^n^± and iVg± (T) = ^ng± (T) and us- 
ing Eq. ( ^ , we rewrite the rate equation given by Eq. ( p3| ) 
in the form 



dN, 



a d 



dt ^ -^^(-^e±^'■)+aF(^^e'±(T)-A^,^±)/7^(36) 



where the equilibrium temperature T has been obtained 
from Eqs.dl) and (|2£ 



5. A simple special relativistic treatment of the expansion 
of the PEM-pulse 

In the following, using the definitions introduced in Sec- 
tion 4, we study the general properties of PEM-pulses. As 
the zeroth order approximation, we do not consider: (i) 
the gravitational interaction since it can be considered a 
perturbation against the total energy of the PEM-pulse 
and we describe the expanding fluid by a special relativis- 
tic set of equations, (ii) the radiation energy lost during 
expansion since this is a very small fraction of total energy 
before the dense PEM-pulse becomes transparent. 

Analogous to Eq.(^), from Eq.(|l|), we have, in the 
general case of presence of baryonic matter, along each 
flow line. 



(29) djuBV) 



dr 



0. 



(37) 



For a small shell of volume AT^ at radius r, for the ex- 
pansion of that shell from the volume AVo to the volume 
AV, we obtain 



AV 



AV7(r) 



AVo AVo7o(r)^ 



(38) 



where AV is the volume of the shell in the coordinate 
frame and it relates to the proper volume AV in the co- 
moving frame by AV^ — 7(r)AV, where 7(7") defined in 
Eq.(^9|) is the 7-factor of the shell at the radius r. 

Similarly from Eq.(|2^), using the equation of state 
(pT]), along the flow lines we obtain 



dlne-l-rdlnt/ = 0. 



(39) 



With a very small variation of the volume from AVo to 
AV and the internal energy density e along the flow lines 
we obtain 



AVY 
AVo 



AV 
AVo 



lo{r) 



(40) 



where the thermal index F given by ( 21 ) is computed in for 
each value of e, p as a function of AV. The thermal index F 
is a slowly- varying function of the state with values around 
4/3. 

We know that as the PEM-pulse expands, the over- 
all energy conservation requires that the change of total 
internal proper energy of a shell in the fluid will be com- 
pensated by a change in the bulk kinetic energy of the 
shell. We clearly have for the total internal proper energy 
of the PEM-pulse in the comoving frame 



£av= / pdV^ pdV, 

JAV JaV 



(41) 



where dV{dV) is the proper (coordinate) volume of a fluid 
element in the shell of the PEM-pulse, and p is defined 
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by (^l]). The change of the kinetic energy of a fluid ele- 
ment in the shell of the PEM-pulse due to an infinitesimal 
expansion is given by, 

dK^[-f{r)^l]{dE + pBdV). (42) 
The total energy conservation for that shell implies 



£aVo - ^Av = / PodV - / pdV = / dK, 
'AVo "'AV JaVo 



(43) 



where £av is the total internal proper energy of the shell 
including baryon masses. 

In order to model the relativistic expansion of plasma 
fluid analytically, we assume that E and D, as defined 
by Eq.(^), are constant in space over the volume AV. 
As a consequence, we can describe the system in terms 
of the relativistic Lorentz factor j{r), the volume AV in 
coordinate frame and the quantities e and pB known in 
the comoving frame. We can rewrite the left-hand-side of 
Eq.(|3|) as 

(eo + p°Bho{r)AVo - (e + PBh{r)AV, (44) 
and the right-hand-side of Eq. ( ^3| ) as 
(e+PB)AV7(r)(7(r)-l)-(eo+P^)AVo7o(r)(7o(r)-l).(45) 
Thus Eq.(||) simpfifies to 

(eo + p?j)7o (OAVo = (e + PBh\r)AV, (46) 
which leads the solution 



7(0 = 7o(0l 



'(eo+p^)AVo 
{e + pB)AV ■ 



(47) 



In the special relativistic case the rate equation (P3| 
for reduces to 



dt 



= -4|-(^'"e±7(0^'') 



(48) 



In the shell of the PEM-pulse, we define the coordinate 
number-density of in equilibrium to be N^± (T) = 
7(r)ne± (T) and correspondingly the coordinate number- 
density of to be Nf,± = j{r)ne± ■ With the assumption 
that ne± in Eq.(|2^) is constant in space over the volume 
AV, we obtain the equation for the evolution of the 
coordinate number-density as they fall out of equilibrium 
with the PEM-pulse, as seen by a coordinate observer 



d ^ .r 1 dAV 1 



{Nli{T)-Nli) .(49) 



Eqs.(^, (^, (^) and (|49|) are a complete set of equa- 
tions describing the relativistic expansion of the shell in 
the PEM-pulse. 

If we now turn from a single shell of the PEM-pulse 
to a finite distribution of shells over the PEM-pulse, we 
introduce the average values of the proper internal-energy. 



baryon-mass, baryon-number and pair-number densities 
{e,pB,nB,ne±), and E = je, D = jpB, N^± = ^{r)n^± 
for the PEM-pulse, where the average relativistic 7-factor 
is defined by, 

l = ^J^lir)dV, (50) 

and V is the total coordinate volume of the PEM-pulse. 
Considering the volume expansion of the PEM-pulse from 
Vo to V, we transform Eqs.(||,|^) into 

V V7 



^ ^ _ 

flB Vo Vo7c 

14 



X r / - \ r 
V \ / 7 



Vo 



7o 



(51) 



(52) 



where V = 7V is the total proper volume of the PEM- 
pulse, following from the definition 7 in Eq.(pO|). Analo- 
gously, Eq.(|4^) becomes 

(eo+P?j)7oVo = (e--fpB)Vf , 



and Eq.(H^) becomes 



7 = 7o 



(eo+P?j)Vo 



Pb)V 



(53) 



(54) 



(55) 



while Eq.(p9[) becomes 

I (7V,± )^-N,^i^+mj^ [Nl, iT) - Nl^ ) , 

where the coordinate number-density of e^-pairs in equi- 
librium is iVe± (T) = 7ne± (T) and the coordinate number- 
density of e^-pairs is iVgi = ^n^± . For an infinitesimal 
expansion of the coordinate volume from Vo to V in the 
coordinate time interval t — to, we can discretize the dif- 
ferential Eq.(|55|) for numerical computations. 



6. Three models of different geometry for the expansion 
of the PEM-pulse 

We analyze three possible patterns of expansion of the 
PEM-pulse for those cases where the relationship between 
the average relativistic 7 factor and the radial component 
of the four- velocity Ur{r) is determined by the following 
respective conditions (see Eqs.( p9|j50| )) 

— Spherical model: assuming the radial component of the 
four- velocity Ur{r) — U^, where U is the radial com- 
ponent of the four- velocity at the moving outer surface 
(r = 7?.(t)) of the PEM-pulse, the 7-factor and the ve- 
locity Vr are 



7 



3 



2u{i + u^)^ - uii + u^y- 



7 



- \n{U + Vl + lP) 



(56) 



this distribution expands keeping an uniform density 
profile, decreasing with time, similar to a portion of a 
Fricdmann Universe. 
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— Slab 1: assuming U{r) — Ur ~ const., the constant 
width of the expanding slab V — Ro in the coordinate 
frame of the PEM-pulse, 7 and Vr are 



7 



Vr = 



Ur_ 

7 



(57) 



this distribution does not need any averaging process. 
Slab 2: assuming a constant width R2 — Ri = Ro of 
the expanding slab in the comoving frame of the PEM- 
pulse, 7 and Vr are 



7= + Vr 



Ur 

— ■) 

7 



(58) 



This distribution needs an averaging procedure and 
i?i < r < R2, i.e. f is an intermediate radius in the 
slab. 

In these three specific models of relativistic expansion, the 
average relativistic gamma factor 7 is related differently 
to the radial component of the four-velocity U{r) and the 
velocity of expansion. As a consequence, this leads to dis- 
tinct relationships between the coordinate radius TZ and 
the coordinate time t in different models. This leads to 
the monotonically increasing 7-factor as a function of the 
coordinate radius (or time) having a distinct slope (see 
Fig.|^). This slope varies from model to model. In princi- 
ple, we could have an infinite number of models by defin- 
ing the geometry of the expanding fluid. Thus, to find 
out which expanding pattern of PEM-pulses is the most 
physically realistic, we need to confront the results of our 
proposed theoretical simplified models with the numerical 
results based on the hydrodynamic Eqs.(^,^,|5|) and see 
if any one of them is in agreement with these results. 

Eqs. ( ^ , ^2|j5^ , |55| ) can be used to make numerical inte- 
grations to relate the quantities 7,6, pb and N^± to 70, eo, 
p°g and N°± in terms of the ratio ^ . We can consider each 
step of the adiabatic process of expansion from volume Vo 
to volume V to be infinitesimal. Given initial values of the 
quantities 70 , £0 , Pg and N°± in the volume Vo and tem- 
perature To as well as N°±{To), we can obtain the values 
of 7, e and pB corresponding to the volume V by using 
computer to solve these three dynamical Eqs.(|l 52 5^). 
With the value of 7 obtained we find the corresponding 
velocity Vr by Eqs.(^,^,|7 58) and the coordinate time 
t = to + At, At = ATZ/Vr, where to,t are the times cor- 
responding to expanding volumes Vo , V respectively and 
TZ — TZo + ATZ, where STZ is the variation of the coordi- 
nate radius at the surface of the PEM-pulse as its volume 
changes from volume Vo to volume V. Then we find the 
corresponding number of pairs Nq± by solving the rate 

and 



N 



equation (p5D. Next we evaluate the fraction ^ ^ 
the appropriate Fermi integrals in the right-hand-side of 
Eqs.( p^ , p5[ ), so that we can determine the temperature 
T of the equilibrium state corresponding to the volume V 
from the proper energy-density through Eqs.(p4,Eq). Then 



by iterative computing steps, Vo,Vi,V2,- • •, we can de- 
scribe the whole process of the expanding PEM-pulse. As 
a result, given initial values 7 — I, TZo ~ rds, To, P°b, Co 
and 7Vg± = N^±{To) of a dyadosphere discussed in sec- 
tion 1, we can obtain 7(7?.) and energy-density e(TZ) of 
the expanding PEM-pulse in terms of the coordinate ra- 
dius TZ and time t. The calculations are continued until 
the plasma becomes optically thin as expressed by the 
condition (p8[). In fact, Eqs.(|53|,p^ are equivalent to 



dt 



= 0, 



(59) 



the conservation of the total energy of the PEM-pulse at 
each expanding step, and we check that it is fulfilled at 
each numerical iteration step. 

From the conservation of entropy it follows that 
asymptotically we have 



{VT^)T<mc- ^ n 
{VT^)T>rnc- 4 



(60) 



exactly for the same reasons and physics scenario dis- 
cussed in the cosmological framework by Weinberg, see 
e.g. Eq. (15.6.37) of Weinberg 1972: the present equation 
follows by substituting the comoving volume factor R^ by 
the comoving volume V in the cited equation. The same 
considerations when repeated for the conservation of the 
total energy p^V = P7^V following from Eq. (^9|) then 
lead to 



{VT^-i) 



T<mc^ 



(VT^j) 



T>mc^ 



11 

T 



(61) 



The ratio of these last two relations then gives asymptot- 
ically 



To = (T7)t>„c2 = (^7) 



T<mc^ J 



(62) 



where To is the initial average temperature of the dyado- 
sphere at rest, given in Fig.|Q. Eq.(|6^) explains the ap- 
proximate constancy of T7 shown in Fig.[0. 

7. Analysis of the Spectrum and Lightcurve at the Time 
of Decoupling 

We are interested in the observed spectrum at the 
time of decoupling. To calculate the spectrum, we as- 
sume that (i) the plasma fluid of coupled e+e" pairs 
and photons undergoes adiabatic expansion and does 
not emit radiative energy before they decouple (this 
assumption is relaxed in Preparata fc Ruffini 199"9| and 



Preparata, Ruffini & Xue 1999); (ii) the e+e^ pairs and 
photons are in equilibrium at the same temperature T 
when they decouple. Thus the photons in the fluid frame 
(denoted with a prime: ') are described by a Planck dis- 
tribution of the form 



,{T') 



exp[ 



(63) 



8 



Remo RufBni et al.: On the Pair Electromagnetic Pulse of a Black Hole with Electromagnetic Structure 



10 ' 



10 ■ 



10 ■ 



^~ 1000 



100 



10 




10 s 



10'" 



10" 
Radius (cm) 



101=^ 



10' = 



Fig. 3. Lorentz gamma factor 7 as a function of radius. 
Three models for the expansion pattern of the PEM-pulse 
are compared with the results of the one dimensional hy- 
drodynamic code for a lOOOAf© black hole with charge to 
mass ratio ^ = 0.1. The 1-D code has an expansion pat- 
tern that strongly resembles that of a shell with constant 
coordinate thickness. 



where T' is the local temperature in the fluid frame of one 
thin shell. However, u^/e^ is a relativistic invariant result- 
ing from the Liouville theorem (see e.g. Ehlcrs 1971 ). This 
implies e/T is also a relativistic invariant. So a Planck 
distribution in an emitter's rest-frame with temperature 
T' will appear Planckian to a moving observer, but with 
boosted temperature T = T'/(7(l — | cos 9)) where ^ cos 9 
is the component of fluid velocity directed toward the ob- 
server. Thus 



Ue{9,V,r) 



exp{'j{l 



1^) 



(64) 



gives the observed spectrum of a blackbody with fluid- 
frame temperature T' moving at velocity v and angle 9 
with respect to the observer. 

We wish to calculate the spectrum from a spherical, 
relativistically expanding shell as seen by a distant ob- 
server. We know the velocity v and proper temperature T' 
and shell radius R, so we integrate over volume and angle 
with respect to the observer. We thus get the observed 
number spectrum N^, per photon energy e, per steradian, 
of a relativistically expanding spherical shell with radius 
i?, thickness dR in cm, velocity v, Lorentz factor 7 and 
fluid-frame temperature T' to be (in photons/eV/47r) 



N,{v,T',R) 



f dV— = (5.23 X lO'^UwR^dR— 
J e 



10 • 



5000 - 



1000 



500 




100 1000 
Mass/Mg 

Fig. 4. In the expansion model of a shell with constant 
coordinate thickness, the decoupling gamma- factor 7gnai 
as a function of EMBH masses is plotted with charge to 
mass ratio = 0.1. 



log 



1 - exp[—fe{l + f )/r] 
1 - exp[—fe{l - f )/T'] 



(65) 



which has a maximum at emax — I.397T' eV for 7 ^ 
1. We may then sum this spectrum over all shells (the 
zones in our computer code) of our PEM-pulse to get the 
total spectrum. Since we had assumed the photons are 
thermal in the comoving frame, our spectrum has now an 
high frequency exponential tail, and the spectrum appear 
as not thermal. In Fig.||] we see the observed number 
spectral peak energy as a function of the EMBH mass 
with charge to mass ratio ^ = 0.1. 

To obtain the observed light curve e{t) we again de- 
compose the spherical PEM-pulse into concentric shells 
and consider two effects. First is the relative arrival time 
of the first light from each shell: light from outer shells 
will be observed before light from inner shells. Second is 
the shape of the light curve from a single shell. 

Emission from the moving PEM-pulse is beamed along 
the direction of travel within an angle 9 ^ I/7. The sur- 
face of simultaneity of a relativistically expanding spher- 
ical shell as seen by an observer is an ellipsoid (see e.g. 
Fenimorc et al. 1996). The time of intersection of an ex- 



panding ellipse and a fixed shell of radius R as a function 
of 9 (i.e. the time at which emission from this intersection 
circle is received) is 



t — —(1 cost 

V c 



R 



(1 - cos9)R+-{l 



V ' c' 



(66) 



Wc find that, integrating our boosted Planck distribution 
of photons Eq.(^4|) over frequency, a relativistically ex- 
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1e+10 




Radius (cm) 



Fig. 5. The number-densities {Ne± /cm~^) of pairs ng± 
(solid line) as obtained from the rate equation (p3|), and 
n^± (T) (dashed Hne) as computed by Fermi-integrals with 
zero chemical potential, provided the temperature T de- 
termined by the equiUbrium condition (|2^), are plotted for 
a IOOOMq black hole with charge to mass ratio ^ = 0.1. 
For T < m^c^, the two curves strongly diverge. 



panding shell of radius R and thickness dR will have a 
time profile (in ergs/second/47r) 



e{T,v,T',R) 



a fT'\^cR^dR 



(67) 



for (1 - f ) < r < (1 + f ), where r = and a is the 
radiation constant. Thus the larger the radius R of the 
expanding shell at the moment of emission, the broader 
the tail on the light curve. The final light curve is 
then constructed by summing the signal from all shells. 

A useful quantity which may be compared with ob- 
servations is tgo; the time over which 90% of the emis- 
sion is received. The tgo duration for a single shell can 
be calculated in the relativistic limit (7 ^ 1) by defining 
e(T) = eoT~'* and setting 



1-7+T90 



1-^ 



s{T)dT ^ 0.9 



1+^ 



e{T)dT, 



(68) 



which gives Tgo 
seconds is 



4.3/7^. So the emission duration in 



^90 



Tgo — = 4.3^- 
C J'^C 



(69) 



for 7 1. Burst durations are shown in Fig.|Q] for a range 
of EMBH masses with charge to mass ratio ^ — 0.1. 



Radius (cm) 



1e+12 



Fig. 6. The reheating ratio T^V/T^Vo defined by Eq.(|o|) 
is plotted as a function of radius for a lOOOM© black hole 
with charge to mass ratio = 0.1. The rate equation (23) 
with definition ( p5| ) naturally leads to the value after 
e~'"e~-annihilation has occurred. 



8. Results and Discussion 

Preliminary results have been presented in Ruffini 
(1999a), Ruffini et al.(1999) and Ruflini (1999b). In the 
results presented here we consider a range of black hole 
masses from 1 to W^Mq and we assume a charge of 
Q = 0-lQmax = 0.1). Over this mass range, and with 
this charge, a large amount of energy (~ 10^^ — 10^'' ergs) 
is available to be deposited within the dyadosphere, in 
the form of e~^e~ pairs, due to the process of vacuum po- 
larization (Fig.[0). In Fig.|^ correspondingly the average 
energy of the created pairs are given as a function of the 
EMBH mass. These energies and spectral distributions are 
in the range of interest for gamma-ray bursts. 

We have analyzed the evolution of this PEM pulse by a 
simple relativistic treatment in flat space for three differ- 
ent cases: 1) an expansion with a spatial component of the 
four velocity proportional to the distance, 2) an expansion 
with a width of the PEM pulse constant in the coordinate 
space and 3) an expansion with a width of the PEM pulse 
constant in the comoving frame. We have then compared 
and contrasted these results with the general relativistic 
hydrodynamic calculations. In Fig.[|j we see a comparison 
of the Lorentz factor of the expanding fluid as a function 
of radius for all of the models. We can see that the one- 
dimensional code matches the expansion pattern of a shell 
of constant coordinate thickness. The astrophysically un- 
precedented large Lorentz factors are reproduced in Fig. [4] 
as a function of the EMBH mass. 
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Fig. 7. Temperature of the plasma is shown over a range 
of masses. Tjnit is the average initial temperature of the 
plasma deposited around the EMBH, Tobs = 1.397gnail" 
is the observed temperature of the plasma at decoupling 
while T' is the comoving temperature at decoupling. No- 
tice that 7T' ~ Tinit as expected from Eq. (|6|) . 
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Fig. 8. The peak of the observed number spectrum as a 
function of the EMBH mass is plotted with charge to mass 
ratio ^ = 0.1. The observed peak temperature Tobs of the 
number spectrum is calculated from Eq. ( |65| ) . The temper- 
ature for the slab of constant coordinate thickness is given 
by Tobs = 1.397fi„aiT'; the peak of Eq.®. 



In Fig.|| and Figs.||,§ we describe the approach to 
transparency, as given by Eq.(p8|). In Fig.|^ we plot the 
number-densities of pairs n^± given by the rate equation 
( p3| ) and rigi (T) computed from a Fermi-integral with zero 
chemical potential with the temperature T determined by 
the equilibrium condition (^J). It clearly indicates that 
the pairs fall out of equilibrium as the temperature 
drops below the threshold of e^-pair annihilation. As a 
consequence of pair e^-annihilation, the crossover of this 
reheating process is shown in Fig. ^ . Fig. illustrates the 
extreme relativistic nature of the PEM pulse expansion: at 
decoupling, the local comoving plasma temperature is < 1 
keV, but is boosted by 7finai to '--^ 1 MeV and its relation 
to the initial energy of the PEM pulse in the dyadosphere 
is presented. 

The spectrum and light-curve of the emitted pulse 
of radiation are analyzed as described in Section ^ We 
see the results for both the hydrodynamic and slab cal- 
culations compared in Figs.||] & The 1-D hydro- 
dynamic code spectral peak and burst duration tgo are 
consistently somewhat lower than those of the slab cal- 
culations due to gravitational and hydrodynamical ef- 
fects. The peak photon energy of the photon number 
spectrum is in the range of 100 — 1000 keV, which is 
consistent with observed gamma-ray bursts. The dura- 
tion of the burst is characterized by the tgo time, the 
time taken to receive 90% of the radiation. Over the 
mass range analyzed, tgo is of order a millisecond. This 
is shorter than observed gamma-ray bursts 1 — 10 




Fig. 9. The duration of the emission at decoupling is rep- 
resented by igo, the time over which 90% of the emission 
is received, plotted over a range of black hole masses. The 
1-D hydrodynamic igo is calculated from Eq. (|67|) and the 
slab fgo is calculated from Eq. (^) . Note that the 1-D and 
slab models deviate for black hole mass M = IMq. This 
is because the energy density Eq.(0) deposited in the 1-D 
computations is not well approximated by a constant over 
the volume as in the slab model for small black hole mass. 
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seconds, see e.g. Tavani 1998), but the collision of the 
PEM pulse with the remnant left over by the process of 
gravitational collapse produces longer emission times, see 



e.g. 


Ruffini et al. 1999b, 


Preparata & Ruffini 1999 


Preparata, Ruffini & Xue 199£. 



We have demonstrated that the concept of an EMBH 
in vacuum can produce a gamma-ray burst having many 
of the general characteristics of observed bursts validating 
the scenario presented in Ruffini 199^ This simplified 
model merits further detailed study of the collapse of 
highly magnetized stars in galaxies (M ~ IOMq) and 
gravitational collapse process occurring in galactic nuclei 
(M ~ 10* — IO^Mq). In particular, it is important to iden- 
tify how, during the process of gravitational collapse to 
an EMBH, the breakdown of the magneto-hydrodynamic 
conditions leads to the formation of the dyadosphere. Fu- 
ture work will explore some additional important aspects 
of the PEM pulse including (i) the presence of baryonic 
matter left in the remnant around the EMBH during 
the process of gravitational collapse and its effect on 
the relativistic expansion (Ruffini et al. 1999b), (ii) the 
radiation emitted prior to decoupling and its effect on the 



dynamics of the PEM pulse (Preparata & Ruffini 1999 



and Preparata, Ruffini & Xue 1999), and (iii) the 
emission in the latest stages of evolution of the 
PEM pulse and its comparison with the observed 



bursts and afterglows in (Preparata & Ruffini 1999 and 



Preparata, Ruffini fc Xue 199£) 
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